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ABSTRACT 

In this paper, we introduce the local Minkowski Functions and apply them as the statistics to study 
the local properties of WMAP Cold Spot (CS) at different scales. We find that these local statistics 
always excess at the scale of 5° comparing with other spots of WMAP data, which clearly presents 
the non-Gaussianity of CS and its characteristic scale. Meanwhile, we find that the cosmic texture can 
excellently explain all the excesses in these statistics, which supports the cosmic texture explanation 
of WMAP CS. 

Subject headings: cosmic microwave background radiation — early universe — methods: data analysis 
— methods: statistical 
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1. INTRODUCTION 

Soon after the release of observations of the NASA 
Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite on the Cosmic Microwave Background (CMB) tem- 
perature and polarization anisotropics, some anomalies 
in C MB field have also been reported (Bennett et al. 
l2Qllf ). Among these, an extremely Cold Spot (CS) cen- 
tered at the position (/ = 209°, h = —57°) and a charac- 
teristic scale about 10° in the Galactic Coordinate was 
detected in the Spherical Mexican Hat Wavelet (SMHW) 
non-Gaussian analysis (Vielva et al. 2004). Comparing 
with the distribution derived from the isotropic and 
Gaussian CMB simulations, due to this CS, the SMHW 
coefficients of W MAP data have an excess of kurtosis 
(jCruz et al.ll2"005h . The non-Gaussian CS has also been 
confirmed by using some other statistics (Cruz et al. 
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■ [2QQ5I : INaselskv et al.l 120101 : 
Vielval [2QloF At the same 
time, some analyses on WMAP CS have also been 
done, such as the non- Gaussian tests for the different 
detectors and different frequency channels o f WMAP 
satellite (IVielva et al.l I2004J : [Cruz et al. 2005), the in- 
vestigation of the NVSS sources (Rudnick et al. 2007; 
iSmith fc Huterer 1 |2010[ ), the survey around the CS 
with MegaCam on th e Canada- France- Hawaii Telescope 
(jGranett et "all 120091 ) the redshift surve y using VIMOS 
on VLT towards CS (jBremer et al.|[2QTQ ^. and the cross- 
correl ation between WM AP and Faraday depth rotation 
map (jHansen et al.ll20 f2). 

Since then, various alternative explanations for the 
CS have been proposed, ^ including the pos sible fore- 
grounds (iCruz et al.' 2006: 'Hanse n et al.ll20"l2h . Sunyaev- 
Zeldovich effect j Cruz et al. 2008), the supervo id in the 
Universe (llnoue fc Silk 2006, 2007; Inoue 2012), and the 
cosmic texture (ICruz et al. 2007b, 2008). Due to the 
fact that, nearly all the explanations of CS are related 
to the local characters of the CMB field, the studies on 
the local properties o f CS are necessary. In the previ- 
ous work (|ZhaZI[20Tl; ). we have studied the local non- 
Gaussian properties of WMAP CS by the local mean 
temperature, variance, skewness and kurtosis. We found 
the excesses of the local variance and skewness in the 
large scales with R > 5°, which implies that WMAP 



CS is prefer a large-scale non-Gaussian structure to a 
combination of some small structures. In this paper, we 
shall focus on the same topic by using the different local 
statistics. We introduce the local Minkowski Functions 
(MFs), and apply them to the WMAP data, in particular 
the WMAP CS. We find that comparing with other spots 
in WMAP data, these MFs around CS have significant 
anomalities at the scale with ~ 5°, which clearly shows 
the non-Gaussian signal of CS and indicates its charac- 
teristic scale. However, in the larger or smaller scales, 
the MFs of CS are all quite normal. Meanwhile, similar 
to (Zh^ [20T2l ). we find that these local non-Gaussianity 
of WMAP CS can be excellently explained by a cosmic 
texture. If subtracting this texture from WMAP data, 
all these anomilies of MFs disappear. So our local anal- 
ysis of the CS in this paper also strongly supports the 
cosmic texture explanation. 

The rest of the paper is organized as follows: In Section 
2, we introduce the WMAP data, which will be used in 
the analysis. In Section 3, we define the local MFs and 
apply them into WMAP data. Section 4 summarizes the 
main results of this paper. 

2. THE WMAP DATA 

In our analysis, we shall use the WMAP data includ- 
ing the ILC7 map and the NILC5 map. The WMAP in- 
strument is composed of 10 DAs spanning fiv e frequen- 
cies from 23 to 94 GHz (jBennett et al.ll2003D . The in- 
ternal linear combination (ILC) method has been used 
by WMAP team to generate the WMAP ILC maps 
fflinshaw et al.ll2007l : IGold et al. 2011). The 7-year ILC 
(written as "ILC7") map is a weighted combination from 
all five original frequency bands, which are smoothed to a 
common resol ution of one d egree. For the 5-year WMAP 
data, in ( Delab rouille et all f2009) the authors have made 
a higher resolution CMB ILC map (written as "NILC5"), 
an implementation of a constrained linear combination of 
the channels with minimum error variance on a frame of 
spherical called needletfl. In this paper, we will con- 
sider both these two maps for the analysis. Note that 
these WMAP data have the same resolution parameter 

^ The similar map for 7-year WMAP data is recently gotten in 
(|Basak Delabrouille 2012). 
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^side = 512, and the corresponding total pixel number 
A^pix = 3145728. 

3. APPLYING LOCAL MINKOWSKI FUNCTIONS 
TO WMAP DATA 

3.1. Minkowski Functions 

MFs characterize the morphological properties of con- 
vex, compact sets in an n-dimensional space. On 
the 2-dimensional spherical surface S^, any morpho- 
logical property can be expanded as a linear combi- 
nation of three MFs, which represent the area, cir- 
cumference and integrated geodesic curv ature of an ex- 
cursion set (jSchmalzing &: Gorskil HoOSh . For a given 
threshold z^, it is convenient to define the excursion 
set Qjy and its boundary dQi, of a smooth scalar field 
u as follows: Qjy = G S^|ix(x) > and dQ^ = 
{x G S^|ii(x) = z^j. Then, the MFs vq^ vi and V2 can 
be written as (jSchmalzing fc Gorskil 1 19981 ). 
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where da and d^ denote the surface element of and 
the line element along dQ^^ respectively. And n is 
the geodesic curvature. Given a pixelated map with 
field u{xi)^ these MFs ca n be numerically calculated by 
the f ormulae ( Schmalzing fc Gorskil [T998c iLim fc SimonI 
[20121 ) 

Vi{iy) = ^Y.Uiy,Xk), (^ = 0,1,2), (2) 
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where 



lo{iy, Xk) :=Q{u - v) 
5{u — v) 



Xi{v, Xk) := ^^^^ ^^ Ui{xk), 
T2{v,Xk) := — U2[Xk), (3) 
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Note that u-i denotes the covariant differentiation of u 
with respect to the coordinate i. The delta function in 
these formulae can be numerically approximated through 
a discretization of threshold space in bins of width Az^ 
by the stepfunction 6n{x) = (Az/)-^[6(x + Aiy/2) - 
Q{x — Au/2)]. The expectation values of three MFs for a 
Gaussian random field are also derived in (Tomita 1986 ), 
which have been explicitly expressed in equations (14) 
and (15) in C Schmalzing fc Gorskil [19981 ). 

These MFs have been applied by cosmologists 
to look from derivations fro m Ga ussianity of the 
pertu rb ations in the CMB 
19981 : ISchmalzing fc Gorskil 
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^999"; 'Eriksen et al.' '2QQ4t iHikage et al.l I2QQ6L I2QQ8I: 
Komatsu et al. 2009; Hikage et al. 2009t iMatsubaral 
2OIO1: iHik age fc Mats ubarai |2012i). In particular, in 
(|Lim fc Simon .20 12.) the authors used the MFs to probe 



the cold/hot disk-like structure in the CMB. However, 
they found these statistics are noise-dominated for the 
WMAP CS. For a single WMAP or PLANCK resolution 
map, the method can only detect the highly prominent 
disk, i.e. extremely cold or hot disk with quite large 
area. These can be easily understood as follows: since 
the MFs are constructed on the global sky, the local 
non- Gaussianity of the cold/hot spot in the relatively 
small scale has been diluted to be too small. 

3.2. Local Minkowski Functions 

To study the local p roperties of WMAP CS, similar 
to the previous works ( Bernui fc Reboucas I l2009l l2010l 
2012; Zhao 2012), in this paper we shall define the local 
MFs as the statistics of CMB field. For a given full- 
sky WMAP data with Nside = 512 (ILC7 or NILC5), 
we smooth them using a Gaussian filter with a smooth- 
ing scale of ^5. Since MFs are sensitive to the smooth- 
ing scale of a density field and thereby we can obtain a 
variety of information from density fields by using dif- 
ferent smoothing levels. In this paper, we focus on both 
ILC fields smoothed by six different smoothing scales 10', 
20', 30', 40', 50', 60'. Then, we can construct the cor- 
responding full-sky maps: Uiixk) and U2{xk) defined in 
Eq. (|3|), where u is the corresponding ILC temperature 
anisotropic map. 

Now, we can define the local MFs. Let Vt{Oj^(j)j]R) 
be a spherical cap with aperture of R degree, centered 
at {Oj,(t)j). The local MFs Vi{Tj) {i = 0,1,2) in this cap 
can be calculated by using Eq. which are denoted 
as vj {v; R) in the rest of this paper. But here the sum 
only counts the pixels inside the cap Vt{Oj^(j)j]R). Note 
that in our calculation, the binning range of threshold v 
is set to be -3.0 to 3.0 with 24 equally spaced bins oiv/a 
{a is the standard deviation of i^-field in this cap) per 
each MF. In order to quantify the same kind of MFs by 
a single quantity, following (Hikage et al.ll2009f ), for each 
i (the type of MFs), {j, R} (the cap) we can define the 
as follows: 

x[viii^^r,R)-vl''{,ya';R)], (4) 

where a and a' denote the binning number of thresh- 
old values. vf^{u; R) is the theoretical value of vj (z/; i?), 
which is independent of the superscript j. S is the corre- 
sponding covariance matrix. Although, in principle the 
theoretical values vf^{i'; R) for the random Gauss ian field 
can be calculated by the analytical formulae (|Tomital 
1986; Schmalzing fc Gorski 1998), there are systemati- 
cal differences from the numerical results due to the bin- 
ning of threshold ( Lim fc Simon 2012). In this paper, 
we avoid this problem by replacing the theoretical value 
vf^{h'; R) by the quantity {vj [v] i?)), which is the average 

value of all vj{u; R) with \bj\ > 30°. And the covariance 
matrix E can also be numerically calculated by these 
quantities vl{v]R). Note that the Galactic plane have 
been excluded to reduce the effect of foreground resid- 
uals. Hereafter, we denote this quantity as X^{R). 

Clearly, the values X^{R) obtained in this way for each 
cap can be viewed as a measure of non- Gaussianity in 






Figure 1. Xo(R) maps (upper), Xi(R) maps (middle), and 
X2(R) maps (lower) for ILC7 data with Os = 60^ smoothing. In 
the left panels, we have used = 2°, and in right panels, = 5° 
is chosen. 



Figure 2. Xo(R) maps (upper), Xi(R) maps (middle), and 
X2(R) maps (lower) for NILC5 data with R = 2° . In the left 
panels, the NILC5 map is smoothed by 6s = 10^, and in right 
panels, Os = 40^ smoothing is applied. 



the direction of the center of cap {Oj^(j)j). For a given 
aperture i?, we scan the celestial sphere with evenly dis- 
tributed spherical caps, and build the Xq{R)-^ Xi{R)-, 
X2(i^)-maps. In our analysis, we have chosen the loca- 
tions of centroids of spots to be the pixels in N^ide = 64 
resolution. By choosing different R values, one can study 
the local properties of CMB field at different scales. 
Same to the y(i?)-map (or S{R)-, K{R)-mdips) defined 

in (jZhaQll2Q12l ). here we also find that (R) always max- 
imize at the edge of the circles, rather than the center 
of circles. To overcome this problem and localize the 
non- Gaussian sources, we define the average quantities 

^/(^):=T^E^'(^)' (^ = 0'1'2), (5) 

where A/'pix is again the pixel number in the j*^ cap. 

We apply the method to the ILC7 data by choosing 
R = 2° and Os = 60^ The Xi maps are presented in 
Fig. [T] (left panels), which clearly show that these local 
statistics, in particular third MF ^2, are very sensitive 
to the foreground residuals and various point sources. 
The Galactic plane is clearly presented, which is the 
non-Gaussian area caused by the foreground residuals 
in ILC7 map. In addition, two important point sources 
at (/ = 209.5°, b = 20.1°) and (/ = 184.9°, b = 5.98°), 
as well as several small ones, also clearly found in the 
Xi maps. So we expect these local statistics with small 
R values can be used to identify the point sources and 
foreground residuals, which will be discussed in a sepa- 
rate paper. If we choose = 5°, from the right panels in 
Fig. [T]we find the similar results, except for some non- 
Gaussianities of small point sources, which have been 
diluted for this larger R case. 

Let us turn to the NILC5 map, which has the much 
higher resolution than ILC7. We firstly study effect of 
different levels of smoothing. By adopting R = 2°, in 
Fig. [2] we plot the Xi maps for Og = 10' (left panels) and 
Os = 40' (right panels). Interesting enough for the lower 



smoothing case, from the left middle and left lower pan- 
els we find the clear structure of Nobs, i-e. the effective 
observations of WMAP for each pixels. We find a smaller 
Nobs, which follows a smaller pixel- noise variance, corre- 
sponds a larger Xf . So, these two local MF statistics can 
also be used to search for the morphology of the pixel- 
noise variance, which has been hidden in the temperature 
anisotropy map. We leave this as a future work. Here, 
in order to reduce the effect of them, we should choose 
a larger smoothing parameter Os as in the right panel 
in Fig. [21 where we find the morphology of the pixel- 
noise variance disappears. But the non-Gaussianity in 
the Galactic place is still there, due to the heavy con- 
tamination caused by foreground residuals. 

Choosing a common smoothing parameter Os = 60', 
in Fig. [3] we plot the Xi maps for R = 2° (left) and 

= 5° (right). By comparing with those in Fig. [H we 
find that although several non- Gaussian point sources 
and the foreground residuals in Galactic plane are still 
there, NILC5 is much cl eaner than ILC7, as claimed in 
(|Delabrouille et al.l[2009h . 

3.3. Local Properties of WMAP Cold Spot 

We now focus on the local properties of WMAP CS 
at (/ = 209°, b = -57°) in Galactic Coordinate. To 
quantify the properties of CS, we define the probabilities 
p- values as follows, 

p,(i?):=Prob(x/(i?)>X,0(i?)), (i = 0,1,2), (6) 

where Xj (R) are the values of spots in the region of |6| > 
30° in Xi{R) map, and Xf{R) is that for CS. Note that 
the Galactic plane has been excluded to remove the effect 
of foreground residuals. So for any given Xi{R) map, a 
smaller Pi{R) value indicates a larger non-Gaussianity 
around CS. 

In Fig. m we plot the probabilities Pi{R) as a function 
of R for ILC7 map with different levels of smoothing. 
We find that Pi{R) > 0.01 is holden for the cases with 





Figure 3. Xo(R) maps (upper), Xi(R) maps (middle), and 
X2(R) maps (lower) for NILC5 data with Os = 60^ smoothing. In 
the left panels, we have used R = 2° , and in right panels, R = 5° 
is chosen. 
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Figure 4. The probabilities Pi(R) versus i?- values for ILC7 based 
Xi maps. From upper left to lower right, the smoothing parameter 
Os = 10^ 20^ 30^ 40^ 50^ eO^ are used. In each panel, black hue 
is for i = 0, blue line is for i = 1 and red line is for i = 2. 



R < 4° and > 6°, which shows that CS in these scales 
is quite normal comparing with other spots in ILC7 data. 
However, when ~ 5°, the Pi{R) values are system- 
atically small for all the three MFs, which means that 
CS causes the significant non-Gaussianity in the scale 
of ~ 10°, which is independent of the smoothing levels 
of ILC7 map. These features can also be found in Fig. 
[TJ where we find CS is quite normal when R = 2° (left 
panels). However, when R = 5° (right panels), are 
significantly large around CS position for all the three 
MF maps. These results are consistent with other au- 
thors' discovery in the SMHW non -Gaussian analysis 
(jVielva et al.l I2QQ4I : iCruz et al.ll2QQ5[ ). So we conclude 
that the local MF statistics can clearly present the non- 
Gaussianity of CS and its characteristic scale. Moreover, 
the normality of CS in the small scales, i.e. < 4°, also 
implies that the WMAP CS is prefer to a large-scale non- 
Gaussian structure, rather than a combination of some 
small structures. Same results are also obtained from 
the NILC5 analysis (see Fig. [5j): so long as 6s > 30' 
is satisfied, where the effect of pixel-noise variances has 
been well suppressed, the CS causes the significant non- 
Gaus sianity at the scale of R ^ 5°. 

In (jCruz et al.ll20Q7bL l2008f ). by studying the temper- 
ature and area of CS, the authors found that the cosmic 
texture, rather than the other explanations, provides an 
excellent interpretation for the WMAP Cold Spot. Here, 
we shall study if the local anomalies of CS found in this 
work is consistent with the cosmic texture explanation. 

The profile for the CMB temperature fiuctuation 
caused by a collapsing cosmic texture is given by 
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(7) 



where is the angle from the center, e is the amplitude 
parameter, and i9c is the scale parameter, i^^ = a/3/2'?^c- 
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Figure 5. Same with Fig. [1 but for NILC5. 
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Figure 6. Same to Fig. [Jl but the cosmic texture structure has 
been subtracted in the original ILC7 map. 
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Figure 7. Same with Fig. [U but for NILC5. 



By the Bayesian analysis, the texture parameters were 
obtained e = T.Stll x 10"^ a nd i9c = ^.dtHdeg at 95% 
confidence (jCruz et al.ll2QQ7bt ). 

In our calculation, we adopt the best-fit texture pa- 
rameters e = 7.3 X 10~^ and i9c = 4.9°, and subtract 
this cosmic texture structure from the ILC7 and NILC5 
maps. Then, we repeat the analysis above by using these 
subtracted ILC maps. The corresponding results are pre- 
sented in Fig. [6] for ILC7 and Fig. [7] for NILC5, where 
we find that comparing with the other spots of WMAP 
data, the CS becomes quite normal, i.e. Pi{R) values are 
alwa ys very cl ose to 1. So, we get the same conclusion as 
in ( Zhaoll2Q12 l: The local analysis of WMAP CS by us- 
ing the local MF statistics strongly supports the cosmic 
texture explanation. 

4. CONCLUSION 

Since the release of WMAP data, many attentions have 
been paid to study the non-Gaussian CS at (/ = 209°, 
b = —57°) in Galastic Coordinate, which might be pro- 
duced by various small-scale contaminations, such as 
point sources or foregrounds, the super void in the Uni- 
verse, or some phase transition in the early Universe. In 
this paper, in order to identify the WMAP CS and dis- 
criminate various explanations, we study the local prop- 
erties of CS by introducing the local MFs as the statis- 
tics. We find that comparing with other spots of WMAP 
data, these local statistics around CS always excess at the 
scale of ~ 5°. So consistent with the discovery of other 
methods, this study clearly shows the non-Gaussianity of 
CS and its characteristic scale. The normality of CS at 
the smaller scales of < 5° also suggests that WMAP 
CS is prefer a large-scale non-Gaussian structure to a 
combination of some small structures. Meanwhile, we 
find that the cosmic texture with a characteristic scale 
about 10°, which is claimed to the most promising expla- 
nation of CS by many authors, can excellently account 
for these anomalies of the local statistics. So our analysis 
supports the cosmic texture explanation for WMAP CS. 
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